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ABSTRACT 

In  this  paper  we  investigate  the  relationship  between  the  stretching 
tensor  J)  and  the  logarithmic  (Hencky)  strain  tnJT,  with  V  the  left 
stretch  tensor,  we  establish  the  simple  formula 

£  -  UnV)°  -  sym  (gnT1)  , 

which  holds  for  arbitrary  three-dimensional  motions.  Here  F  is  the 
deformation  gradient,  (Anjj)0  is  the  time  derivative  of  InV  measured  in  a 
coordinate  system  which  rotates  with  the  left  principal  strain  axes,  and  ft 
is  the  spin  of  the  right  principal  strain  axes.  We  use  this  formula  to  show 
that  £  ■  (In vf,  (or,  equivalently,  Q  *  ,  the  Jaumann  derivative  of 

taV)  ,  if  and  only  if  the  characteristic  spaces  of  the  right  stretch  tensor 
JJ  are  constant  on  any  time  interval  in  which  the  number  of  distinct 
principal  stretches  is  constant.  Finally,  we  show  that  the  asymptotic 
approximation 

£  -  +  0(e3) 

holds  whenever  the  displacement  gradient  g  satisfies  g,g  ■  0(e). 
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SIGNIFICANCE  AND  EXPLANATION 


\ 

Experimentalists  often  report  strain  data  using  logarithmic  measures  of 
strain,  and  these  measures  have  also  been  utilized  extensively  for  theoretical 
purposes  in  the  recent  engineering  literature.  with  V  »  X^f  4)  f ^ 

(summation  implied)  the  left  stretch  tensor,,  we  define  the  logarithmic 
(Hencky)  strain  tensor  by  AnV  -  (AnX.)  f,  ®  f..  The  X.  and  f,  are 

**  X  '''l  ^1  1 

referred  to  as  the  principal  stretches  and  left  principal  axes  of  strain, 
respectively.  'This  logarithmic  strain  measure  has  been  shown  to  have  certain 
advantages  in  the  formulation  of  the  deformation  theory  of  plasticity. 

Whenever  the  principal  axes  of  strain  are  fixed  throughout  the  motion,  it 
is  known  that  the  simple  relationship 

E  -  (An X  )* 

holds,  with  E  the  stretching  tensor.  Previous  work,  however,  has  failed  to 
find  an  analogous  formula  in  the  general  case.  Some  researchers  have 
therefore  concluded  that  the  Hencky  strain  measure  has  very  limited 
applicability . 

The  purpose  of  this  paper  is  to  extend  the  above-,  formula  to  arbitrary 
three-dimensional  deformations^  We  are  able  to  find  a  simple  relationship 
between  E  and  (AnV)° ,  where  fyie  time  derivative  of  Anv  is  now  measured 
by  an  observer  rotating  with  the  r^ .  By  generalizing  the  idea  of  fixed 
principal  axes,  we  show  precisely  when  the  exact  formula  D  «  (AnV)°  holds, 
we  conclude  by  showing  that,  for  sufficiently  small  deformations,  (AnV)*  is 
an  excellent  approximation  to  E* 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


CM  TOE  RELATIONSHIP  BETWEEN  TOE  LOCARITHGC 
STRAIN  RATE  AND  TOE  STRETCHING  TENSOR 

Marten  E.  Gurtin  and  Kathleen  Spear* 

1.  Introduction. 

The  logarithmic  strain  introduced  by  Kencky  1 3]  has  long 
enjoyed  favored  treatment  in  the  metallurgical  and  materials 
science  literature,  where  it  is  referred  to  as  the  "true"  or 
"natural”  strain.  Its  use,  however,  has  been  primarily  limited 
to  studies  in  which  the  principal  axes  of  strain  are  fixed.  In 
such  problems,  the  simple  relationship 


a  -  itojrr 

exists  between  the  stretching  tensor  a  an*  logarithm  of 
the  left  stretch  tensor  £. 

In  this  paper,  we  investigate  the  question  of  whether  an 
analogous  relationship  exists  for  general  three-dimensional 
deformations.  Truesdsll  and  Toupin  [9]*  note  that  the  Hencky 
strain  has  never  been  successfully  applied  in  general.  Recent 
attempts  to  remedy  this  situation  include  the  work  of  Hill  [4] , 
who  argues  that  logarithmic  strain  measures  have  inherent 
advantages  in  certain  constitutive  inequalities.  He  finds  the 
rather  complicated  relation 

P>ii'  i  - 

(4n  U)*^  -  L2(XiXj4n(li/kj)Dij)/uJ-k2)f  i  #  j, 

*pp.  269-270. 
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with  the  components  of  (to  y)  *  and  £  taken  with  respect 
to  eigenvector  bases  of  y  and  y,  respectively.  (Here  U,  is 
the  right  stretch  tensor  and  the  are  the  principal  stretches; 

i.e.,  the  eigenvalues  of  U.) 

Both  Hutchinson  and  Neale  [S]  and  Storen  and  Rice  [8]  find 
Hencky  strain  to  be  useful  in  the  formulation  of  the  deformation 
theory  of  plasticity ,  although  Storen  and  Rice  decide  that  the 
general  relation  between  D  and  (toy)*  is  "very  complicated". 
They  conclude  that  toy  is  "essentially  intractable”  as  a  strain 
measure.  Fitzgerald  [1]  decides  that  the  utility  of  logarithmic 
strain  is  limited  to  problems  with  fixed  principal  strain  axes. 

In  addition  to  this  negative  concensus  on  the  applicability  of 
Hencky  strain  measures,  the  above  authors  all  consider  only 
strains  for  which  the  principal  stretches  are  distinct,  and 
do  not  rigorously  treat  the  possibility  that  the  principal  axes 
of  strain  may  not  be  uniquely  defined. 

We  here  attempt  to  give  a  complete  answer  to  the  general 
three-dimensional  problem.  We  begin  by  establishing  a  simple 
general  formula  relating  y  and  to  V: 

a  “  (to  xr  -  ■yn>(ysry"1) . 

Here  £  is  the  deformation  gradient,  is  the  spin  of  the  right 
principal  axes  of  strain,  and  (toy)*  is  the  time  derivative  of 
to  y  measured  by  an  observer  rotating  with  the  left  principal 
axes  of  strain.  We  use  this  formula  to  prove  that  y  ■  (to  V)° 

(or,  equivalently,  y  ■  (toy)*,  the  Jaumann  derivative  of  toy), 
if  and  only  if  the  characteristic  spaces  of  y  are  independent 
of  time  on  any  time  interval  in  which  the  number  of  distinct 


principal  stratchaa  is  constant.  Finally,  we  show  that  for 
motions  which  ara  small  in  tha  sansa  that  tha  displacement  gradient 
£  ■  £  -  J,  satisif as  *  0  ( €) ,  ^  •  0  ( t) ,  we  have  tha  estimate 

Ji  -  (*n  jy*  ■»*  0 ( €3) . 

Since  the  error  term  is  two  ordars-of-magnitude  higher  than 
(in  V)*,  wa  conclude  that  (in  jy*  is,  in  fact,  an  excellent 
approximation  to  Q,  in  this  instance. 


2.  Preliminaries. 

Consider  a  motion  of  a  body,  and  let  F(t)  denote  the 
deformation  gradient  corresponding  to  a  given  material  point. 

At  each  t  (in  a  fixed  time  interval)  F(t)  is  a  3x3  matrix 
with  strictly-positive  determinant,  and  hence  admits  the  polar 
decomposition 

i  -  sa  -  a. 

where  y,(t)  and  £(t)  are  symmetric,  positive-definite,  while 
&(t)  is  proper  orthogonal.  £  and  £,  respectively,  are  called 
the  right  and  left  stretch  tensors.  We  assume  that  £  is  smooth 
(i.e.,  continuously  differentiable) j  then 


k-  faf1 


is  the  velocity  gradient. 


£  -  sym  J.  -  fife**?) 
the  stretching  tensor,  and 

3(t-kT) 

m 

the  spin  tensor.  (Here  denotes  the  transpose  of  £.) 

The  principal  stretches  are  the  eigenvalues  X^,  i  •  1,2,3. 

2 

of  a  (or  yj ;  since  £  smooth,  we  may,  without  loss  in 
generality,  choose  the  three  functions  X,  (t)  to  be  smooth3 


e  follow  the  notation  and  terminology  of  [2] 
Cf.  [2],  p.  23. 

Cf.  Kato  (6],  Thm.  6.8,  p.  111. 
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in  t.  An  orthonormal  basis  1^3  of  eigenvectors  of  U,  (where 
corresponds  to  X.)  is  called  a  right  principal  basis «  We  assume 
that  one  such  basis  is  given,  and  that  each  r^(t)  is  smooth 

in  t.  Then  [j^3  defined  by 


is  a  left  principal  basis;  that  is,  (4^3  ia  an  orthonormal  basis  of 
eigenvectors  of  £  corresponding  to  (X^l.  In  view  of  the  spectral 
theorem,  we  have  the  representations 

H.  *  *i£i  ®  Si*  X  “  *i£y.  •  ii* 

Here  £  ®  £  with  components  a^b^  is  the  tensor  product  of  a 
and  b,  and  summation  over  repeated  indices  is  implied.  Also, 

£  *  xi£i  •  £±* 

The  characteristic  space  for  U  at  time  t,  corresponding  to 
the  principal  stretch  X^,  say,  is  the  set  of  all  vectors  £  such 


JL(t)X-  XA(t)v, 

so  that  coincident  principal  stretches  correspond  to  the  same 
characteristic  space. 

Let  tj^3  be  a  smooth,  time-dependent  orthonormal  basis.  Then 
the  corresponding  twirl  tensor  Q  is  the  skew  tensor  function 
defined  by 

^This  is  an  assumption:  £  smooth  does  not  necessarily  yield  the 
existence  of  a  smooth  basis  lj^3  (cf.  Kato  [6],  Example  5.9,  p.  115) 

2 

Cf e  9  €«9*  (  [2]  ,  P«  11* 


ii  “  £6d- 

(0  is  skew  since  Si’Sj  “  6ij  i®pli«*  that  « 

S^*  Sj  +  Sj.*  Sj  “  ® • )  Given  any  smooth  tensor  function  A,  the 
co-rotatlonal  derivative  A*  of  A  relative  to  £e^}  is  defined 
by 

ar  -  i  *  as  -  & 


and  represents  the  time  derivative  of  A  measured  by  an  observer 

rotating  with  (t) ] .  Another  important  notion  is  the  Jaumann 
1  * 

derivative  A  of  A,  given  by 

£  -  k  ♦  &  - 

with  W  the  spin. 

Finally,  the  tensor  logarithm,  <tn,  maps  symmetric,  positive- 
definite  matrices  into  symmetric  matrices  and  is  defined  to  be 
the  inverse  of  the  exponential  function.  In  particular. 


*n  to  ty£i  ®  ri.  In  V-  ttn  X.)^  ® 


T 


Cf.,  e.g.,  [7],  p.  155 
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3.  Relationship  between  D  and  liny). 

ffi^ojce^n  J,.  Let  ^  denote  the  twirl  tensor  corresponding 
to  Cj^J  and  (  )•  the  co-rotational  derivative  relative  to 
C^J.  Then 

B  -  (in  -  sym<| IQ,.?1) .  (4) 

Proof.  Differentiation  of  <3)2  gives 

u»xr  -  cVVfc  ®4l  +  (to  Vi  ♦  (to  V&  ®  i 
-  ( vvfe  •  h  *  St(to  »  -  . 

with  the  twirl  tensor  corresponding  to  Thus, 

<4«  X)*  - 

Next,  by  (1)  and  (2) , 

fc  -  <Vi  ®  £i  +  xii  ®  £i  +  Xiii  ®  i> Ej  ® 

-  <  Wii  *  k  + 

-  <*»  ♦,&  -  ttyf1*  (5) 

This  completes  the  proof,  as  (4)  is  the  symmetric  part  of  (5). 

Since 

(4n  -  (4n  X?  *  +  <4n  X)W  -  W(*n  y) , 

(4)  may  also  be  written  in  the  form 


Z  tp 


D  ■  (to  V)  +  (to  V) (n^-W)  -  (n^-W) (to  V)  -  syrnfFn^)  ;  (6) 
we  will  use  (6)  in  later  calculations. 

Remark.  The  term  FQ^F"1  represents  the  spin  of  the  right 
principal  basis  (r^}  as  measured  by  am  observer  deforming  with 
the  body.  (Let  tj^}  be  a  basis  fixed  in  space,  and  a  basis 

deforming  with  the  body,  i.e.,  h.  “  fv^.  Then  the  components  of 
relative  to  (j^]  are  the  same  as  those  of  FQ^f"*  relative 

Remark.  Similar  arguments  show  that 

a  -  £(to  Ji)‘ST  -  »ym(£SrE-1), 

where  (  )*  is  the  co-rotational  derivative  relative  to  [r^], 
rather  tham  [j^J. 
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4.  When  does  D  -  Un  V)°? 

In  this  section  we  present  a  condition  on  the  principal  strain 
axes  under  which  the  formula  D  ■  (in  jy*  is  valid.  To  state  this 
result  precisely  it  is  necessary  to  extend  the  notion  of  fixed 
principal  axes  of  strain.  This  idea  makes  no  sense  when 
two  or  more  of  the  principal  strains  coalesce,  for  at  those 
times  the  axes  are  not.  uniquely  defined.  The  characteristic 
spaces,  however,  are  uniquely  defined,  but  change  in  type  depending 
on  the  number  n(t)  of  distinct  principal  stretches:  when 
n(t)  -3  the  characteristic  spaces  are  three  mutually  perpen¬ 
dicular  lines;  when  n(t)  >2  the  characteristic  spaces  are  a 
line  and  a  plane  perpendicular  to  it;  when  n(t)  >1  the  char¬ 
acteristic  space  is  all  of  JR"*  .  Thus  it  only  makes  sense  to 
demand  that  the  characteristic  spaces  be  independent  of  time 
on  time  intervals  during  which  n(t)  remains  constant. 

Precisely  then,  let  n(t)  denote  the  number  of  distinct 
rlnclpal  stretches  at  time  t.  We  say  that  the  right  principal 
axes  are  essentially  fixed  if  the  right  stretch  tensor  U  has 


characteristic  spaces  which  are  fixed  in  time  on  any  time 
interval  during  which  n(t)  is  constant. 

We  are  now  in  a  position  to  state  our  main  result. 


yftfpjyw  2.  The  following  are  equivalent: 

(a)  The  right  principal  axes  are  essentially  fixed. 

(b)  The  co-rotational  derivative  of  In  V  correspondln 
to  the  left  principal  basis  U. J  satisfies 


(c)  The  Jaumann  derivative  of  to  V  satisfies 


D  -  (to  V)\ 

This  theorem  has  the  following  obvious 

Corollary.  Suppose  that  the  three  principal  stretches  are 
distinct.  Then  the  formulae 

D  ■  (In  V)#  ,  D.(toV)* 

hold  if  and  only  if  the  three  right  principal  axes  are  fixed  for 
all  time. 

Of  course,  the  right  principal  axes  are  the  three  lines 
generated  by  the  basis  vectors  . 

The  next  two  lemmas  facilitate  the  proof  of  Theorem  2.  In 
these  lemmas  and  in  their  proof,  0^,  and  0^,  respectively,  denote 
the  twirl  tensors  corresponding  to  the  right  and  left  principal 
bases,  and 

Lemma  JL.  Let  T  be  a  time  interval  of  nonzero  length  on 
which  n(t)  is  constant.  Then  the  following  are  egui valent: 

(i)  During  T  the  characteristic  spaces  of  are 
Independent  of  time. 


on  T 
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Proof.  We  begin  by  noting  that 

£i  •<«&-£,£>  £j  "  2£i*&£j  ’  £i*SrS&j 

■  (*>i“*j)  Ei*  figXj. 

Cage  1  (n-3) .  By  (7) ,  (ii)  ia  equivalent  to  »  £  on  T, 
and  since  (i)  and  (ii)  are  equivalent. 

Case  2  (n-2) .  Here  (using  the  spectral  theorem)  we  may, 
without  loss  in  generality,  write  £  on  T  in  the  form 

X(t)  -  ®£l(t)  +  S(t)II-Ei(t)«Ei(t)], 

so  that  the  line  spanned  by  £l(t)  and  the  plane  perpendicular 
to  £l(t)  are  the  characteristic  spaces.  Assume  that  (ii)  holds. 
Then,  since  1^  i<  1^,  we  conclude  from  (7)  that 

~jr~l  *  Ea’StSi  "  0  (8) 

on  T.  Thus 

£i*£2  "  £i*£3  '  °» 

and  since  r^jr^  ■  0,  we  have  ■  constant  on  T,  which  implies  (i) . 

Conversely,  if  (i)  holds,  the  above  argument  in  reverse 
shows  that  (9)  is  valid.  Condition  (ii)  then  follows  from  (7) , 
as  X2  “  *3  on  T* 

Case  3  (n«l) .  Here 


a<t)  -  Mt)*, 


I 


on  T,  and  the  conditions  (i)  and  (ii)  are  satisfied  identically. 

This  completes  the  proof  of  Lemma  1. 

£•  The  following  are  equivalent ; 

(i)  The  right  principal  axes  are  essentially  fixed. 

(ii)  («or  all  time) . 

Proof.  That  (ii)  implies  (i)  follows  trivially  from  Lemma  1. 
If  (i)  holds,  then  0^  commutes  with  U  on  all  time  intervals 
of  nonzero  length  during  which  n(t)  is  constant  (by  Lemma  1) . 

In  view  of  the  continuity  of  the  X^(t),  n(t)  is  piecewise 
constant,  and  for  any  time  tg  at  which  n(t)  jumps  there  are 
right  and  left  intervals  (a,tg)  and  (tg,b)  of  nonzero  length 
on  which  n(t)  is  constant.  Continuity  of  and  a  then 

gives  (ii) . 

Proof  of  Theorem  2.  (a)  •  (b) .  By  (4) ,  D  •  (4n  £)•  is 

equivalent  to  svm(Ffltfl~1)  -  0.  Since 

iSirl1  *  EfS^RT^E** 

the  latter  condition  is  equivalent  to  the  requirement  that 

W1  -  a"1fir5i  -  o, 

or  equivalently 

a2Sr " 


(9) 


As  is  known,1  a  symmetric  tensor  A  commutes  with  a  tensor  B  if 

and  only  if  g,  leaves  invariant  the  characteristic  spaces  of  A. 

2 

Thus,  since  the  characteristic  spaces  of  £  and  a  coincide, 

(9)  is  equivalent  to  (ii)  and  hence  (i)  of  Lemma  2. 

(a)  .  •  (c) .  Assume  I>  ■  (to  Then,  by  (6), 

<*n  V)  (j^-W)  -  (&-W)  (to  V)  -  symt^jf1)  -  0. 

Since  W  -  -  BkwtgJ^g,-1) , 

(to  VJskwJgC^F-1}  -  skw^f1)  (to  X>  -  symtF^E,"1)  -  g,. 

Using  V  -  RjjR"  and  the  isotropy  of  to 

R[  (to  y)skw(ynru~1)  -  skwtyn.u"1)  (to  jy  -  syraju^ji-1)  ]g,T  -  o, 

or  equivalently, 

(to  a>  (sisru’1+ji'*1firJi)  -  <jtwr1*ar1&0 (to  "  W1 

Since  U  commutes  with  to  a» 

a)  (u^+sirU2)  -  <a2jat+&a2> (Ul  a>  “  a2flc  -  qji2-  <10> 

As  before,  we  will  show  that  U  commutes  with  by  showing 

that  leaves  invariant  the  characteristic  spaces  of  a*  Thus 

fix  the  time  t,  let  \  denote  a  principal  stretch  and  A  the 
corresponding  characteristic  space  for  a  ”  £(*) *  choose  jr  €  A, 
and  let 

&  *  flrtt* 


Applying  (10)  to  a  and  noting  that 


JJ,  Jt  -  X  yt,  (to  ft)*  -  (to  X)j», 


we  arrive  at 

um  k  -  c*  u#  (jiVx2x>  -  ji2^  -  x2& .  (u 

Let  3  be  a  principal  stretch  with  3  #  X  and  let  e  be  a 
corresponding  eigenvector,  so  that  £-w  *■  0.  Then  taking  the 
inner  product  of  (11)  with  e  yields 

(to  3  -  to  X)  (32+X2)  (2.-3E)  -  (32-X2)  (a.jt). 

Thus  either  -  0  or 


and,  as  we  shall  show  at  the  end  of  the  proof,  (12)  implies  that 
3  •  X.  Thus  £•&  *  0.  We  have  shown  that  x  is  orthogonal  to 
all  characteristic  spaces  except  A.  Hence  x  €  A  and 
leaves  invariant  the  characteristic  spaces  of  ft. 

Conversely,  if  commutes  with  ft,  then  0^  commutes  with 

to  ft.  Since  £  commutes  with  to  ft,  equation  (10)  holds  trivially, 
and  reversing  the  arguments  leading  to  (10)  yields  the  condition 
ft-  (to  ft)\ 

Thus  (c)  is  equivalent  to  (ii)  and  hence  (i)  of  Lemma  2. 

To  complete  the  proof  we  have  only  to  show  that  (12)  implies 
X  -  0.  Thus  let  p  -  to(3/X).  Then  (12)  is  equivalent  the  equation 
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2 

Since  the  derivative  of  tanh  p  is  sech  p  £  1,  with  equality 
holding  only  at  p  *  0,  the  graphs  of  p  and  tanh  p  intersect 
only  at  p  »  0.  This  completes  the  proof  of  Theorem  2. 

Remark.  Some  other  relations  which  hold  when  the  right 
essentially  fixed  are 

<<"&)•-  Bf1  -  rt  -  &TR£- 


principal  ax< 
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as  I A  |  -•  0.  In  (16)x,  A  is  arbitrary;  in  (X6)  2,3'  h  ia 
symmetric. 

Our  next  step  will  be  to  estimate  the  right  aide  of  the 
identity 

(to  V)*  -  (to  jy  *  (to  jyw  -  £  to  (17) 

Let 


£  -  sym  ii  ■  il™ 

Then  by  (13) , 

x2  -  sx,T  - i *  m +  bb,t 

and  (16) 2t2  Yiald 

X  -  i  ♦  X  ♦  i<SST-£2>  ♦ 

X  •  X  ♦  T  »*  -  £2  ♦  £<£> »  £<a>  -  o<€3)  • 


(18) 


To  derive  an  asymptotic  expansion  for  (to  jy  *  w«  write  £(£) 
for  to  v  with  a  considered  as  a  function  of  £.  Then  (18)  2 
is  the  Taylor  expansion 

IW  -  £<jy  ♦  £'<jy  W  ♦  §  £"(£>  (£*£]  +  C(£),  (19) 


where  £'  (£)  and  £"  (£) .  are  the  first  and  second  (Frechet) 
derivatives  of  a  at  2*  with  £'  (0)  (£1  linear  in  £* 
t"(0)  I&x'Ba1  *ynMtric  and  bilinear  in  <ai'£2>*  zf  wa 
differentiate  (19)  with  respect  to  t,  we  arrive  at 

£'(£)(£]  -  £' (£)[£!  +  £"(£)(£,£)  ♦  £(£>*  • 


On  the  other  hand,  if  we  expand  X'  (H)  about  g,  ■  0 ,  using  (14) 

we  get 


v  <&>  ih]  -  <a>  chi  +  r  (<y  ia,&]  +  o<  €3> . 

Bence 

ccai-  -  o(€3). 


and  differentiating  (18) 2  with  respect  to  tine  yields,  after  some 
work, 

(4a  V)  *  -  £  -  sym(gtl)  +  G£  -  +  0(€3) .  (20) 

Next, 

fc«  of1  - 

•  k -  fi£  +  otc3) * 


and  so 


•  •  3 

£  ■  sym  Jj  *  E  -  sym(gH)  +  0 ( €  ) , 
N  -  skw  Jt  ■  £  +  0(62). 


(21) 


The  estimates  (17),  (18) 2,  (20),  and  (21)  imply  the  desired  result  (15) 


Remark.  Since  £  ■  0(€)  and  (4n  £)*  -  0(€),  the  asymptotic 

*  3  * 

expansion  £  ■  (tn  J£)  +  0(€  )  shows  (tn]()  to  be  an  excellent 

e 

approximation  to  D  when  &  and  &  are  small. 
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